Influence of Carbon Concentration on the Superconductivity in MgC^jNis 
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The influence of carbon concentration on the superconductivity ( SC ) in MgC^Nis has been 
investigated by measuring the low temperature speciflc heat combined with flrst principles electronic 
structure calculation. It is found that the specific heat coefficient 7„ = Can/T of the superconducting 
sample {x ~ f) in normal state is twice that of the non- superconducting one [x ~ 0.85). The 
comparison of measured 7„ and the calculated electronic density of states ( DOS ) shows that the 
effective mass renormalization changes remarkably as the carbon concentration changes. The large 
mass renormalization for the superconducting sample and the low Tc( 7K ) indicate that more than 
one kind of boson mediated electron-electron interactions exist in MgCxNia. 

PACS numbers: 74.25.Bt, 74.20.Rp,74.72. Jt 



The recently discovered superconductor MgCNia I"] 
with an anti-perovskite structure has attracted renewed 
attention, because it is a new intermetalic compound 
following the remarkable discovery of MgB2 and is 
thought to be very close to ferromagnetism (FM) due 
to a large proportion of Ni in each unit cell 3] . Exper- 
iments have shown that the superconductivity (SC) in 
MgC2;Ni3 is very sensitive to the content of carbon and 
disappears below about x = 0.88 Q. Moreover, Ni-site 
doping with Cu (electron doping) and Co (hole doping) 
have very different effect . 

Many efforts have been focused on the origin of SC 
in MgCNis. Some researchers explored the possibihty 
of conventional phonon-coupled pairing 6,7]. However, 
various experiments and theoretical calculations have not 
given a consistent coupling constant A in the frame of 
electron-phonon coupling (EPC) mechanism. For exam- 
ple, A K, 1 — 1.6 was estimated from the calculated plasma 
frequency {up) and measured resistivit y lal , while the 
measured upper critical field Hc2 0, liol lllj suggested a 
much larger value of A « 2.5 — 3. Furthermore, two dis- 
crepant values of A = 0.7 1] and 2.2 ^3 were determined 
in the similar specific heat measurements by taking dif- 
ferent approximations. Besides these confusing results, 
it has been proposed that EPC is possibly not the only 
contribution to the coupling mechanism in MgCNis sys- 
tem [§, 10]. Opposed to the conventional pairing mecha- 
nism, it was suggested that FM and SC may coexist due 
to the ferromagnetic instability caused by the high den- 
sit y of electronic states (DOS) at the Fermi energy (Ei?) 
[Sl, [l^ • NMR experiment confirmed the presence of sub- 
stantial ferromagnetic spin fiuctuations in MgCNia [l3| . 
while it also exhibited a clear coherence peak of the nu- 
clear spin-lattice relaxation rate just below Tc, which is 
typical for an isotropic s-wave superconductor. More re- 
cently, Voelker et al. 0] proposed that MgCNia could be 
a multiband sup erconductor with a conventional phonon 
mechanism Il5^ similar to MgB2 0]. Thus, the origin 



of SC in MgCNia is still controversial. One of the most 
illuminative methods is to study the infiucnce of element- 
substitution on SC. 

In this Letter, we compare the specific heat data of 
a superconducting MgCajNia with x « 1 and that of 
a non-superconducting one with x « 0.85. Using the 
tight-binding linear muffin-tin orbital (TB-LMTO) band 
method and coherent potential approximation (CPA) 
[T? . 18], we investigate the carbon concentration depen- 
dence of the electronic structure of MgCajNia . Combining 
the results of the specific heat and band calculations, we 
propose a self-consistent picture to understand the re- 
duced SC in non-stoichiometric MgC^jNia. 

Poly-crystalline samples of MgCa;Ni3 were prepared 
by powder metallurgy method. Details of the prepara- 
tion were published previously 0. The superconduct- 
ing transition of the stoichiometric sample( SCI ) with 
Tc — IK is shown in Fig. ^ A non-superconducting 
sample (NSC08) was synthesized by continually reduc- 
ing the carbon component until the diamagnetism was 
completely suppressed. X-ray diffraction (XRD) patterns 
show nearly single phase in both samples. The lattice pa- 
rameters a determined from XRD are 3.812A and 3.790A 
for SCI and NSC08, respectively. By comparing these 
lattice parameters with the reported a ^ x relation 0], 
we could estimate the carbon content as a; = 0.977 and 
0.850 for SCI and NSC08, respectively. According to 
the Tc ~ a; relation presented in Ref 4], x = 0.977 cor- 
responds to Tc ~ and x = 0.850 is very close to the 
critical stoichiometry of SC in MgCa:Ni3, which is in good 
agreement with our experimental results. The tempera- 
ture dependence of the normalized resistivity is shown in 
Fig. n The heat capacity data presented here were taken 
with the relaxation method |2| based on an Oxford cryo- 
genic system Maglab, in which the heat capacity is deter- 
mined by a direct measurement of the thermal time con- 
stant, T = {Cs + Cadd)/ i^wi here Cg and Cadd are the heat 
capacity of the sample and addenda, respectively, while 
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FIG. 1: The temperature-dependent resistivity (normalized 
by the value at 300K) of the samples SCI and NSC08. Inset: 
The enlarged view of the superconducting transition in AC 
susceptibility measurement. 

the thermal conductance between the chip and a 
thermal link. During the measurement the sample was 
cooled to the lowest temperature under a magnetic field 
( field-cooling ) followed by data acquisition in the warm- 
ing up process. The addenda includes a small sapphire 
chip, a tiny Cernox temperature sensor, small amount 
of Wakefield thermal conducting grease and gold leads 
((/)25/im). All contributions to the heat capacity from the 
addenda have been measured separately and subtracted 
from the total specific heat before further analysis. 

The typical low temperature specific heat at various 
magnetic fields up to 12 Tesla is presented in Fig. 13 For 
sample NSC08, all the data fall into an universal curve, 
indicating the independence of its specific heat on the 
magnetic field. For sample SCI, all the data above Tc are 
independent on the field, and the normal state extends 
to the whole temperature region investigated here when 
H > lOT. The upper critical field Hc2{T) can be deter- 
mined from the position of the specific-heat jump, and 
Hc2{0) ~ IIT is then derived by extrapolating Hc2{T) to 
OK. In order to describe the specific heat data of NSC08 
and the normal state data of SCI, we use the following 
expression: 

CniT) = -fnT + ^T^ + ST^ (1) 

where the linear-T term is the electronic contribution 
with 7„ as Sommerfeld parameter, the second term rep- 
resents the phonon contribution according to the Debye 
approximation, the last term is required to include devi- 
ations from the linear dispersion of the acoustic modes in 
extended temperature range. The fitting results are pre- 
sented in Fig. [51 and the obtained parameters are shown 
in Table ^ Besides the obvious difference between the 
Debye temperatures (6_d) of SCI and NSC08 (due to 
the larger elastic modulus of NSC08 than that of SCI), 
a more prominent difference is found between the Som- 



FIG. 2: A plot of C/T vs for MgCNig (SCI) and 
MgCo.ssNia (NSC08) in various magnetic fields. The solid 
lines are fits to Eq. ^ as discussed in the text. The spe- 
cific heat jump at the superconducting transition tempera- 
ture is completely depressed by the field above lOT. Upper 
inset: The low-temperature specific heat for SCI at various 
magnetic fields from OT to 12T. Lower inset: Fitting for the 
zero-field specific heat data (open circles). 



TABLE I: Fits to Eq. for samples SCI and NSC08. 
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merfeld parameters 7„ of these two samples, this is one 
of the main results in this work. 

Before further analysis, we at first investigate the spe- 
cific heat data of the superconducting sample SCI in zero 
field, as shown in the lower inset of Fig. 12 The solid curve 
represents the theoretical fits, the part above is the fit- 
ting result to the normal state data of SCI as mentioned 
above, while the part below T^, is an attempt to describe 
the data using the specific heat expression of conven- 
tional superconductor: Ces = Aexp{-A{0)/kBT)/T^-^, 
in which ^4 is a function of zero-temperature energy gap 
A(0) and the normal-state DOS g{0)- This equation is 
a low temperature approximation and is not adequate 
above 0.8Tc where A(r) has substantially deviated from 
A(0). This conventional expression can fit our data be- 
low O.STc « 5K quite well, though the returned A(0) 
(1.46meV) leads to a large value of 2A{0) /kgTc « 5 com- 
paring with the weak coupling value 3.5. This indicates 
that MgCNia is a strong coupling superconductor, being 
consistent with the large value of AC/jnTc ~ 1.7 ( the 
expected value for weak coupling is 1.43 ), where AC is 
the specific heat jump at Tc. 

As shown above, the specific heat coefficient 7„ of SCI 
is almost twice that of NSC08. In the frame of strong 
coupling mechanism, given the value of 7„ and N(Ep), 
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an effective mass renormalization fact or ( A) can be esti- 
mated from the following relation [21I 

7n = Y^l{l + X)N{EF) (2) 

In order to obtain N{Ep), we have calculated the elec- 
tronic structures using a TB-LMTO method based on 
Green function formalism within the atomic sphere 
approximation (ASA). The parameterized exchange- 
correlation potential of Ref [2^ and the scalar-relativistic 
Dirac equation l23l| were employed. The coherent poten- 
tial approximation(CPA) used in our calculation allowed 
us to study the electronic structure of MgC^Na for any 
Carbon concentration. In the crystal structure of cu- 
bic antiperovskite-type MgCNa, the atoms occupy the 
positions Mg (0,0,0), C (0.5,0.5,0.5), and Ni (0.5,0.5,0), 
(0.5,0,0.5), (0,0.5,0.5). Self-consistency was reached us- 
ing 1000 k-points within the irreducible wedge of the sim- 
ple cubic Brillouin zone. In the calculations, we adopted 
the experimental determined relation between the lat- 
tice parameter a and the carbon concentration x. A fixed 
radii proportion as 2.078 : 1 : 1.617 (Mg:C:Ni) was taken 
for the stoichiometric situation (x = 1), and only the 
radii of carbon-site was reduced with decreasing x. 

The DOS of MgC^^Ns for x = 1 are presented in the 
upper inset of Fig. El which is consistent with the reports 
of other groups 0, M EE IH ■ The tt* antibonding states 
of Ni-Sd and C-2p are located just below Ep, yielding 
a high DOS peak whose height is proved to be sensitive 
to the exact method. Fortunately, the obtained value of 
N{Ep) = 4.53 eV~^cell~^ is not far away from the re- 
sult of general potential linearized augmented planewave 
method (4.99 eV-^ceir^) |i] and that of full-potential 
LMTO method (4.57 eV-icell-i) Hi, which is a good 
starting point for the following consideration of the dop- 
ing effect. 

Fig. shows the evolution of the DOS around Fermi 
level with the change of carbon concentration x in 
MgCa;Ni3. It is found that the Fermi level slightly shifts 
towards low energy with the decreasing x, while the Ni3d- 
C2p antibonding peak is depressed remarkably. In the 
nearest-neighbor approximation, the degeneration of this 
DOS peak indicate that the carbon-site vacancies locally 
break the hybridization between Ni-3(i and C-2p orbitals 
and hence lead to a redistribution of the electronic states. 
Atom-resolved DOS shows that the near-Ep DOS peak 
of Ni-3(i component and that of G-2p component do not 
depart from each other with decreasing x, indicating 
that carbon-site vacancies does not destroy the Ni3(i- 
C2p bond for all carbon concentrations investigated here. 
Therefore, the decrease of the peak possibly means the 
enhanced itinerancy of electrons. As shown in Fig.O the 
depressed DOS peak just below Ep results in a notable 
reduction of N{Ep). The dependence of N{Ep) on the 
carbon concentration is presented in Fig. 01 We can see a 
obvious decrease of N{Ep) from 4.26eV~^ at a; = 0.977 
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FIG. 3: Electronic DOS of MgCajNis for various carbon con- 
centration calculated within the ASA. The dotted line indi- 
cate the Fermi level. Upper inset: Electronic DOS of MgCNis 
in a lager energy scale. Lower inset: Redistribution of DOS 
with the change of carbon concentration. The thick solid line 
and thin one corresponds to a; = 0.977 and 0.85, respectively. 
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FIG. 4: The dependence of N{Ef) of MgC^Nia on carbon 
concentration. N{Ef) for fixed a = 3.812A is also calculated 
as a comparison. The solid lines are guides to eyes. 

to 3.14eV^^ at x = 0.85. Moreover, the doping depen- 
dence of N{Ep) is linear above a critical doping Xc « 0.88 
which is very close to the carbon content in NSC08 and 
also the lower limit of the carbon concentration for SC 
in MgCajNia 0. Below Xc, the doping dependence of 
N(Ep) becomes much weaker than that above Xc- 

By inserting the experimental value of 7„ and the cal- 
culated N{Ep) into Eq(|21), we can determine the cou- 
pling constant as A = 2.37 and 1.35 for the samples SCI 
and NSC08, respectively. The former is close to that re- 
ported by Walte et al. [l^, while the remarkable differ- 
ence in the coupling strength between a superconducting 
sample and a carbon-deficient non-superconducting one 
is for the first time found in the MgC^jNis system. How- 
ever, applying A = 2.37 to the McMillan's formula [l^ 
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Tc — Y^^~^ gives rise to Tc > SOK, which seems to 
be too large as compared to experimental value Tc = IK. 
In other words, Tc = 7K corresponds to A = 0.42, 
which is much smaller than the value determined here. 
This contradiction can be resolved if we adopt the pic- 
ture that there are two kinds of boson-mediated electron- 
electron interactions, for example, electron-phonon cou- 
pling (EPC) and spin fluctuation (SF) Due to the 
different effect of the spin fluctuations on mass renormal- 
ization and superconducting properties, it should add in 
the mass renormalization term and subtract in the pair- 
ing term for SC. So in the frame of strong coupling, 
the specific heat determined A has two contributions: 
A = Xph + Aspin- Similarly, the simplified McMillan's 
Tc formula should be modified as Isl 



level with the decreasing carbon concentration. Further 
analysis indicate that there is more than one kind of 
boson-mediated electron-electron interactions existing in 
MgC^Nis system. In the frame of strong coupling theory, 
a substantial depression of the electron-phonon coupling 
caused by the decrease of carbon concentration is for the 
first time found in this system and seems to be responsi- 
ble for the impairment of its superconductivity. 
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dation of China (NSFC 19825111, 10074078), the Min- 
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nkbrsf-gl999064302 ), and Chinese Academy of Sciences 
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Using the determined parameters for sample SCI 
(namely, A = 2.2, Tc = 7K and 6 = 351if), one 



obtains A- 



■ph 



A 



spin 



= 0.95 from Eq.( EJ. Combining 



this result with Xph + As 



2.37, we determine the 



strength of both EPC and electron-paramagnon coupling 
as Xph = 1.66 and Xspin = 0.71. Firstly, by assuming that 
Xspin is invariable with the decreasing carbon concen- 
tration, we could obtain Xph = 0.64 for sample NSC08. 
This suggests a strong dependence of the EPC strength 
on the carbon concentration in MgC^rNis. For NSC08, 
Xph = 0.64 has become smaller than Xspin — 0.71, which 
is closely associated with the disappearance of SC. Sec- 
ondly, on the assumption that x = 0.85 is the criti- 
cal stoichiometry for SC, i.e., Xph = Xspin, we obtain 
Xph = Xspin = 0.665. This indicates a 6% decrease of 
Xspin from X = 0.977 to a; = 0.85. However, Xph is lost 
near 60% in the same system, which seems to be respon- 
sible for the suppression of SC in MgC^jNia. 

We have been aware of some recent calculations [l^E^ 
suggesting that both Co doping and Cu doping in 
MgCNis lead to a reduction of N{Ep), which seems to 
be critical in depressing SC in this system (MgCNis^^Yx 
[Y=Co or Cu]). However, these investigations were based 
on the ordered substitutions, i.e., x=l, 2 or 3, which is 
somewhat away from the key problem because the solu- 
bility of Cu is limited to 3% {x = 0.1) in technology and 
the bulk SC disappears for only 1% Co doping {x = 0.03). 
Our further calculations did not show obvious reduction 
of N(Ef) in this doping range. Therefore, some neces- 
sary investigations on the variation of the EPC in this 
system must be supplemented. 

In summary, we have investigated the specific heat 
data in MgC^^Nis system and found a remarkable dif- 
ference between the Sommerfeld parameters of a super- 
conducting sample (x w 1) and a non-superconducting 
one {x « 0.85). Band calculations reveal a distinct de- 
crease of the density of the electronic states at Fermi 
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